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The elimination of decoherence of two-state quantum systems interacting with a thermal reservoir 
through an external controllable driving field is discussed in the present paper. The restriction 
equation with which the external controllable driving field should agree will be derived. Based 
on this, we obtain the time-development equation of the off-diagonal elements of density operator 
in the supersymmetric multiphoton two-state quantum systems, which is helpful for studying the 
polarization evolution in this two-state quantum model. 
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I. INTRODUCTION 

Recently, an area called quantum computation, which involves computers that use the ability of quantum systems 
to be in a superposition of many states, attracts extensive attention of many researchers. However, it is not yet clear 
whether quantum computers are feasible to build [1]. One reason that quantum computers will be difficult, if not 
impossible, to build is decoherence. In the process of decoherence, some qubit or qubits of the computation become 
entangled with the environment, thus in effect "collapsing" the state of the quantum computer [1]. In literature, 
there may exist three schemes to reduce the decoherence: (i) error-avoiding codes [2]; (ii) error-correcting codes 
[3]; (iii) decoherence-avoiding scheme [4]. The third approach to the suppression of decoherence can be realized 
by eliminating the interaction between the (two-state) quantum system and environment (such as a noise field, 
bath, thermal reservoir and so on) in the presence of an external controllable driving field. In this paper, we will 
study the maintenance of coherence via the decoherence-avoiding scheme, and the time evolution of polarization in a 
supersymmetric multiphoton two-level model. On considering the latter problem, we assume that the environmental 
effect on the quantum system under consideration has been eliminated by using the decoherence-avoiding scheme, 
which enables us to treat the polarization evolution problem in the multiphoton model more conveniently [i.e., the 
polarization evolution problem will be investigated under the assumption that the decoherence of the quantum systems 
in the noise field has been reduced). 

In section II, we introduce Zhang's method for treating the decoherence problem [5]. In section III, by employing 
Zhang's approach to the multiphoton two-state quantum system we obtain the time-development equation of the 
off-diagonal element of density operator in this two-state system. 

II. MAINTENANCE OF COHERENCE 

In this section, we review one of the formulation for dealing with the decoherence-avoiding scheme, which was 
suggested by Zhang [5] more recently. The reason for the detailed reappearance of Zhang [5] below is as follows: (i) 
first and foremost, Zhang's formulation is helpful for treating the polarization evolution problem in the multiphoton 
model; (ii) in this paper, the elimination of decoherence is a prerequisite for simplifying the polarization evolution 
problem in the multiphoton model, namely, in studying the polarization evolution problem in the multiphoton two- 
state system, we assume that the entanglement of the quantum system with the environment {e.g., thermal reservoir) 
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has been eliminated by an external driving field. So, we need not consider the decoherence problem of the multiphoton 
system in section III. 

Let us first consider the following model, the Hamiltonian of which takes the form (in the unit U = 1) 



H(t) = y<t 2 + Wfcoj^Ofc + ^2g k <r z (a k + a{) - -[E(t)a+ + E*{t)a_], 



(2.1) 



which can describe the interaction of a two-level atom with a noise field (thermal reservoir). In this Hamiltonian, luq, 
u>h, &z, 9k and E(t) denote the atomic transition frequency, photon frequency with k- mode, the third-component 
Pauli matrix, the coupling coefficient (of atoms to the thermal reservoir) and the external driving field, respectively. 
By using the unitary transformation 



V(t) = exp 



(2.2) 



one can arrive at the Hamiltonian 



Hi(t) =^2gk<J z dk exp(-ia; fe i) + a\ exp(iwfei) - - [E(t) cxp(iLu t)a + + E*(t) exp(-iw t)a-] (2.3) 
k 

in the interaction picture, where use is made of H\{t) = V^(t) [H(t) — i-§i\ V(t). The density operator of the two-level 
atomic system agrees with 



i^ = [H l (t),p I (t)]. 

Let p q i(t) denote the atomic reducible density operator. It follows that the reducible density operator equals 



Pql (t) = Tr r p!(t) = - f^H^t), W), p I (t')}]dt' . 
Jo 



(2.4) 



(2.5) 



If we assume that the thermal reservoir is rather large, then it can be concluded that the reservoir may not change 
much during the time evolution process of the atom-reservoir system. Thus we have pi(t) ~ Pqi{t) Pri(0) , where p r i(0) — 
cxp(— /3i? )/Trexp(— (3H ) with (3 — l/k B T, H = ^}2 k ^ka\ak- If we take pi(t') ~ p\(t) (Markoff approximation), 
then Eq.(2.5) can be rewritten as 



Pgi(t) 



where the integrand can be rewritten as 



Tr r [i? I (t),[H I (t'),P 9 i(t)pn(0)]]dt', 



(2.6) 



Tr r [Hi(i), [#i(i'),p g i(*)/>n(0)]] - Tr r [i/i(i)ffi(i')P 9 i(i)Pri(0) - #i(*)p g i(i)p r i(0)#i(*') 

- Hi (t')p q i(t) PlI (0)Hi (t) + Pq i{t)pA (O)fl-i (t')Hi(t)]. 



(2.7) 



By the aid of the relations < a\a k >r = rik, < a k a\ >r= rife + 1, n k — l/[exp(hu}k/kQT) — 1], one can arrive at 
through lengthy calculation [5] 



/M</) = -2 [~ 
with A = Ai + A 2 + A 3 + A 4 , where 



Pq i(t)\E(t)\ 2 f exp[iw (* - t')]dt' - p g i(t)A + a z p ql (t)<T z A 
Jo 



(2.8 



a 2/- , , n 1 - exp(-iwfct) v-^ 2- l-exp(wj fe i) 

* 7/. Jr. * 



^3 = £flfc(Sfc + l) 



lUJk 

1 - cxp(ia; fc t) 



—lUJk 

1 - exp(-zcj fc t) 



(2.9) 



2 



Ignoring the Lamb-shift term (i.e., taking the real parts of poi = (0|/6 9 i|l), poi = (0|p 9 i|l), pw = (1 |0) , pw = 
(l\p q i\Q)), we can obtain the time-development equation of the off-diagonal elements of density operator, i.e., 



An 



Pw = -2 



sin(uj t) d 2 
loq 4 

sin(u;oi) d 2 

L0 4 



\E(t)\ 2 +J2^t(2n k + l) 

k 

■\E{t)\ 2 + ^9l^n k + l) 



sin(a> fc t) 



sin(wfci) 



Poi, 



PlO- 



Thus it is readily verified that if the envelope of the external driving field satisfies the following condition 

-nM^ |jB(t)|2+ ^ 2 ^ (2 _ fe + 1) sin(^) 
k 



wo 



0. 



OJk 



(2.10) 



(2.H) 



then one have poi(t) = poi(0) and pw(t) — Pio(0), which means the suppression (elimination) of dccoherence in this 
two-level quantum system interacting with the environment (thermal reservoir) through an external controllable field 
E(t)- 

It should be noted again that the above theory was proposed by Zhang [5] . In the next section we will investigate 
the polarization evolution problem in the two-level supcrsymmctric multiphoton Jayncs-Cummings model by making 
use of Zhang's formulation. 



III. POLARIZATION EVOLUTION IN THE MULTIPHOTON TWO-LEVEL QUANTUM SYSTEM 



The multiphoton two-level system that we will consider in this section is the supersymmetric multiphoton Jaynes- 
Cummings model [6,7], the Hamiltonian of which under the rotating wave approximation is given by 



H = ^y<r z + una) a + g{a)) k a_ + g*a k <r + , 



(3-1) 



where and a are the creation and annihilation operators for the electromagnetic field, and obey the commutation 
relation [a, cr] = 1; a± and a z denote the two-level atom operators which satisfy the commutation relation [<j z , a±\ = 
±2ct± ; g(t) and g*(t) are the coupling coefficients and k is the photon number in each atom transition process; u>o(t) 
and u)(t) are respectively the transition frequency and the mode frequency. 

The supersymmetric structure can be found in this multiphoton two-level quantum model by defining the following 
supersymmetric transformation generators [8,9]: 



N = a f a+ ■ 
Q = (a)) k <r 



k-l 1 

<Jr H = 

2 z 2 





(at) fc o 







o 



a a i = 



N = 



a k {a^) k 
{a)) k a k 



a k 




(3.2) 



It is easily verified that (N, N' ,Q, Q^) form supersymmetric generators and have supersymmetric Lie algebra proper- 
ties, i.e., 



[JV,Qt]=-Qt, {Q\Q}=N, {Q,a z } = {Q\a z }=0 
[Q, u z \ = 2Q, [Q\a z ] = -2Q\ (Qt - Q) 2 = - N ' , 



= 0, [N,Q] = Q, 



(3-3) 



where {} denotes the anticommuting bracket. 

Now let us obtain the Hamiltonian of the above multiphoton Jaynes-Cummings model in the interaction picture by 
using the following unitary transformation 



V(t) = cxp 



- — a z +Lua'a) t 
i \ 2 / 



(3.4) 



and the result is 
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Hi(t) = gexp(-i6t)Q + g*exp(iSt)Q^ (3.5) 
with 5 = kuo — uo . Based on Eq.(2.7), by complicated calculation, one can arrive at 

TMtfi(*)fli(*>9i(*)M0)] - < #i(t)fli(*>«i(*) >r 

= gg* {exp[iS{t' - t)] < QQ t p q i{t) > R + exp[-iS{t' - t)] < Q^Q Pq i{t) > R } , (3.6) 

where Q 2 = (Q^) 2 = is applied to the calculation, and 

Tr^^i^/^O)^')] = < H^p^H^t') > R 

= g 2 exp[-iS(t + f)} < Q PqI (t)Q > R +(g*) 2 exp[i6(t + t')] < Qt PgI (t)Qt > R 

+ [exp[iS(*' - i)] < Qp^Q 1 > R + exp[-^(t' - t)] < Q^p qI (t)Q > R ] , (3.7) 

Tr r [tfi(i')p ? iWPri(0)#i(*)] = 9 2 expH^(t + *')] < Qp fl i(t)<9 >R +(3*) 2 cxp[ i( 5(t + t')] < QWWQ* >R 

+ gg* [exp[-iS{t' - t)] < Q Pq i{t)Q^ > R + exp[i<S(i' - t)} < Q^ Pq i(t)Q > R ] , (3.8) 

and 

Tr r [p 9l (i)pri(0)Hi(t')^i(0] = 99* {exp[-iS{t' - t)} < Pq i{t)QQ^ > R + exp^i' - t)\ < Pq i(t)Q^Q > R } . (3.9) 
Thus it follows from (2.7) and the above four expressions (3.6)-(3.9) that 

Tr r [#i(i), [Hi(t'),p 9 i(t)Pri(o)]] = r 1 + r 2 + r 3 , (3.10) 

where 

Ti = 99* {exp[iS(t' - t)] < {a)) k a k > R cr_cr + + exp[-iJ(t' - t)] < a k {a)) k > R a+a^} p ql (t) 

+ gg*p q i(t){exp[-i5(t' -t)\ < {a)) k a k > R a_a+ + cxp[i5(t' - t)} <a k {a)) k > R( 7 + a_}, (3.11) 

T 2 = -2{.g 2 cxpH£(t + 0] < Q Pq i(t)Q > R +(. 9 *) 2 exp[^(i + t')] < QW(*)Q f > R } (3.12) 

and 

T3 = -99* {exp[iS(t' - t)\ + exp[-iS{t' - *)]} (< Q Pq i(t)Q^ > R + < Q*p q i(t)Q > R ) . (3.13) 

With the help of the following relations a fc ( a t) fe |m) = |m), < a fe (at) fe > R = ^±j^, < (at)fc a fe > R = 

< (a t ) fc (a t ) fc > R = 0, (7+|+) = 0, cr_ | — ) = 0, (+|<7_ = 0, (-\a+ = 0, one can arrive at 



(0|T 1 |l)= 55 *exp[i<y(t / -t)] 
<0|T 2 |1)=0, 



(m + fc)! ml 



ml (to — fc)! 



<o|p,i(*)|i>, 



(to + fc)! 



(0|T 3 |1) = -gg* {cxp[t5(t> - t)] + exp[-iS(t' - t)}} v ' >' (l|p g i(t)|0). (3.14) 



ml 



If we set p i(t) = (0|p,i(t)|l), pio(t) = (l|p g i(*)|0), poi(i) = (0|p g i(t)|l), and use the equation 

m(t) = - [ '& r <o|[fl- I (t),[/ri(t'),p,i(t)pri(o)]]|i>dt' (3.15) 

Jo 

in accordance with (2.6), then we will get 

Poi(i) = ci(t)poi(t) - c 2 (i)pi (i) (3.16) 

with 



ci(t) = -gg* 



(m + k)l ml 



ml (to — fc)! 



1 — exp(— iSt) , . ^(m + fc)! 

^ i , c 2 (t) = -33* i -L 

io ml 



1 — cxp(— iSt) 1 — cxp(iSt) 
iS —id 



(3.17) 



In the meanwhile, we can obtain the complex conjugation to Eq.(3.16), i.e., 

Pio(t) - c*(t) Pl0 (t) - c 2 (t)p 01 (t). (3.18) 

Thus we obtain the time-development equation of the off-diagonal elements of density operator in this supersymmetric 
multiphoton two-state quantum systems. 

In order to indicate the physical meanings of Eqs.(3.16) and (3.18), we will set u = poi + Pio, v = i(poi — Pio)- 
Eqs.(3.16) and (3.18) can therefore be rewritten as follows 

— = [Reci(i) - c 2 (t)] u + lmci(t)v, 
at 

dv 

— = [Reci(i) + c 2 (t)] v - Imci(i)tt. (3.19) 

It is easily seen that here the physical meanings of u and v are as follows: u and v have close relation to the dispersion 
and the dissipation (absorption/amplification) of polarization, respectively. 



IV. BRIEF DISCUSSION 



To close this section, we briefly discuss the time-evolution equation (3.19) of u and v. Let us consider a simple case 
where the detuning frequency 8 from the atomic transition frequency is vanishing. It follows from the expressions 

(3.17) that d(t) = -gg* (m +, fc)! + (m m! fc) , j t and c 2 (t) = -2gg* {m ^ ] ' t. Note that here a(t) is a real function. So, 

according to Eq.(3.19), we obtain 



i(t) = uq exp 



v(t) = vq exp 



(4.1) 



where A„ = -gg* 



(m-fe)! 



(rn+k)\ 



and A„ = -gg* 



(m-fc)! 



+ 3 



(m+fc)! 



Note that since X u > for k > 0, the 



function u(t) associated with the dispersion will exponentially increase in the process of evolution, while v(t) related 
to the absorption/ amplification will decrease rapidly in this process. In a decoupling case where k = 0, A„ = 0. In 
this case, the real part of polarization (i.e., u) will not vary with the development of the quantum system. However, 
the imaginary part (v) of polarization does not vanish. 
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